Microtubule dynamics depart from wormlike chain model 
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Thermal shape fluctuations of grafted microtubules were studied using high resolution particle 
tracking of attached fluorescent beads. First mode relaxation times were extracted from the mean 
square displacement in the transverse coordinate. For microtubules shorter than ~10^im, the re- 
laxation times were found to follow an L 2 dependence instead of L 4 as expected from the standard 
wormlike chain model. This length dependence is shown to result from a complex length depen- 
dence of the bending stiffness which can be understood as a result of the molecular architecture of 
microtubules. For microtubules shorter than ~5 /im, high drag coefficients indicate contributions 
from internal friction to the fluctuation dynamics. 

PACS numbers: 87.16Ka, 87.15La 
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Microtubules are cytoskeletal protein filaments that 
play an essential role in a multitude of cell functions in 
all eucaryotes. While specialized structures such as cilia, 
flagella and axons require microtubule lengths of up to 
several hundred micrometers, fundamental tasks like cell 
division and intracellular transport involve microtubules 
with lengths that are comparable to or smaller than typ- 
ical eucaryotic cell sizes (10-20 /Ltm). In solution, micro- 
tubules undergo strong thermal shape fluctuations but, 
due to their rigidity, maintain an average direction and 
are therefore referred to as semiflexible. The standard 
model for semiflexible polymers is the wormlike chain 
[H, 0] which envisions a homogeneous, isotropic, continu- 
ously flexible rod characterized by its bending stiffness, k. 

Recently, evidence has accumulated suggesting that 
the wormlike chain as the standard model for semiflexible 
polymers may constitute an oversimplification in the case 
of microtubules due to their highly anisotropic molecu- 
lar architecture 0, 0| • Microtubules consist of strings of 
a- and /3-tubulin heterodimers, so-called protofilaments, 
that are arranged in parallel forming a hollow tube of 
25 nm external diameter. Neighboring protofilaments are 
shifted relative to each other, giving rise to helicity. Both 
the rise per turn and the number of protofilaments may 
vary, however 3 monomers per turn and 13 constituent 
filaments are most common in vivo [5j, [6( . This protofila- 
ment architecture makes microtubules a model system for 
the generalized theory of wormlike bundles 0, Q , which 
in addition to various other biological examples also de- 
scribes the mechanics of carbon nanotube bundles. 

Few studies have addressed the dynamics of ther- 
mal shape fluctuations of microtubules even though the 
timescales of these fluctuations affect the timescales on 
which biological functions occur. Caspi et al. |9|] ex- 
amined the transverse mean square displacement of mi- 



crotubules in networks, finding a power law behavior 
Janson and Dogterom [lfj extracted autocorre- 
lation times for several modes from a shape analysis 
of growing grafted microtubules, while Brangwynne et 
al. [ll| applied this procedure to stabilized fluorescent 
microtubules. For technical reasons, all of these studies 
restricted themselves to lengths of several tens of microm- 
eters, while microtubule lengths relevant to cell functions 
such as cell division are much shorter. In a study of equi- 
librium position distributions jij, the most significant 
systematic deviations from the wormlike chain model 
were found for microtubules shorter than ~20 ^m which 
coincides with the length regime most crucial to cell di- 
vision and intracellular transport. This suggests that, in 
this length regime, corresponding deviations should exist 
for dynamical parameters such as the relaxation time. 

In this letter we present an analysis of first mode relax- 
ation times for microtubules of length 2-30 fim. Micro- 
tubules are grafted to a substrate at one end and data 
is extracted from the mean square displacement of the 
transverse coordinate. Capturing dynamics for micro- 
tubules that are only several /im long is challenging since 
the corresponding spatial and temporal scales become in- 
creasingly small. We overcome these difficulties by using 
small fluorescent beads attached to the microtubules as 
tracer particles. Their position can be tracked at frame 
rates of up to 30 Hz using a high quantum efficiency CCD 
camera. In contrast to shape analysis techniques that 
provide low resolution (tens of nanometers T3, HI) P°~ 
sition data for the whole filament, the use of the tracer 
bead yields spatial information with a precision of a few 
nanometers for one specific point on the filament's con- 
tour. Standard semiflexible polymer models are then 
used to infer relaxation times, bending stiffnesses and 
drag coefficients. 
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In the wormlike chain model, relaxation times for a 
filament are given by [l2} 



(L 4 



(1) 



where £, L, and n denote the drag coefficient per unit 
length, the filament's contour length and its bending stiff- 
ness, respectively, while q n is a mode-dependent numer- 
ical factor which corresponds to a wave number normal- 
ized by length and is determined by the boundary condi- 
tions. An estimate of £ can be obtained from the formula 
for a rigid cylinder of length L and diameter d held steady 
in a homogeneous flow in a liquid of viscosity 77 [l3T ] : 
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hx(L/d) + 2 ln(2) - 1/2' 
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The Stokes friction coefficient for the attached bead is 
7bcad = Qt^W ~ 1-9 x 10~ 9 Ns/m for r\ f» 10 -3 Ns/m and 
r = 100 nm being the viscosity of water and the radius of 
the bead, respectively. The drag Co-^ acting on the mi- 
crotubule is significantly larger than 7b oa d and only be- 
comes comparable for microtubules shorter than ~1 /im 
2l|. Hence, in our case, the attached bead does not need 
to be considered in the analysis of the microtubule dy- 
namics. 

Due to confinement of the filament's position, the mean 
square displacement of the transverse coordinate satu- 
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MSD(f) «V^fl- e-'M U(LJL), (3) 

n 11 

where L a refers to the contour length from the micro- 
tubule attachment point to the bead and f n (L a /L) ac- 
counts for a modulation of the amplitude due to the po- 
sition on the contour. In the case of one free and one 
clamped end, q n = {11 — l/2)n, and the sum is taken over 
all modes. As ((72/91) 4 ~ 80, it follows that the ampli- 
tude of the first mode dominates the mean square dis- 
placement. Furthermore, Eq. |T]) shows that relaxation 
times decrease rapidly with mode number. For fluores- 
cence microscopy studies that rely on integration times 
on the order of several to tens of milliseconds this implies 
that higher modes are generally averaged over during 
imaging and therefore do not contribute significantly to 
the measured MSD. Given the stiffness of microtubules, 
we can therefore safely assume that our measurements 
essentially reflect the behavior of the first mode. Con- 
tributions from the second or higher modes may occur 
during measurements on long microtubules where relax- 
ation times are significantly longer, however these are a 
negligible source of error compared to the statistical un- 
certainties in this length regime. 

Microtubules were polymerized from unlabeled, 
rhodamine-labeled and biotinilated tubulin (Cytoskele- 
ton Inc, Denver, CO) and stabilized in 20 fiM taxol 




FIG. 1: Schematic of the assay. One end of the microtubule 
is covalently grafted on the gold substrate. Lengths are mea- 
sured along the contour starting at the substrate, with L de- 
noting the free end and L a < L marking the position of the 
tracer bead. Different fluorescent labels for the microtubule 
and the bead allow for separate imaging of the two and pre- 
vent photobleaching. 



(T1912, Sigma, St. Louis, MO) to ensure stability of the 
length. Using standard thiol chemistry (see Ref. [4[), mi- 
crotubules were covalently grafted to a microstructured 
gold substrate at one end, leaving the other end free to 
fluctuate in three dimensions. A yellow-green fluorescent 
bead of 200 nm in diameter was tightly bound to the 
microtubule using specific binding between avidin and 
biotin. Fig. [1] shows a schematic of the assembly. For 
19 microtubules of lengths L = 2.2 to 27.9 ^m, 10 000 to 
15 000 frames of the bead's motion were taken at frame 
rates of 10 to 30 Hz using integration times of 18 to 80 ms. 
Two-dimensional position data were obtained from the 
time series using a custom-written Matlab (The Math- 
Works, Natick, MA) single particle tracking algorithm 
involving Gaussian fits to the bead's image. 

MSD data obtained from time averages of finite time 
traces are subject to statistical errors that are nontrivial 
to detect and quantify. In order to avoid biased results 
arising from correlations in the data, error estimates on 
MSD data were obtained as outlined in Ref. (l6j . 

In addition to statistical errors, instrumentation lim- 
itations need to be considered. Finite integration times 
average position fluctuations and thereby result in effec- 
tive low-pass filtering [ljj. We account for the low-pass 
effect by extracting relaxation times from a fit to a cor- 
rected expression for the measured MSD of first mode 
fluctuations [22I ]. 

In Fig. [2] the resulting relaxation time data is plotted 
versus the contour length. For short filaments, the data 
approach a power law r oc L 2 which differs considerably 
from the r a I 4 dependence predicted by the wormlike 
chain model. This observation leads to the conclusion 
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FIG. 2: Relaxation times extracted from MSD. Lines of slope 
4 (solid) and 2 (dashed) are shown for comparison. 

that either the bending stiffness n or the drag coefficient 
£ in Eq. ([1]) must have a hidden dependence on L. 

The persistence length l p — ^(fesT) can be deter- 
mined from the transverse position variance V using 

am 



where a low-pass filtering corrected expression is applied 
[l7j |. Once l p and hence k is known, the drag coefficient 
( can be computed using Eq. (Q]). 

The plots of l p and £ versus L in Figs. and 0] im- 
mediately reveal that both quantities exhibit unexpected 
behavior. The persistence length shows a significant in- 
crease with overall filament length, similar to the length- 
dependence obtained in Ref. [H where, based on the Tim- 
oshenko beam formalism [l^ | that describes macroscopic 
elastic beams with shear contributions, the expression 




was derived. The parameter A accomodates for decreases 
in microtubule stiffness due to shear contributions from 
inter-protofilament sliding. In our case however, a fit 
of Eq. ([5]) cannot provide much information as our data 
range does not include the saturation at predicted for 
very long microtubules. While, for microtubules longer 
than ~5 /im, qualitative agreement of the data with the 
quadratically increasing regime of Eq. §5§ is evident, 
shorter microtubules seem to deviate from this behav- 
ior and level to a plateau. Averaging the values for the 
three shortest microtubules yields a value of 580±100 /im. 
Interestingly, the existence of such a plateau value is pre- 
dicted by the more general theory of wormlike bundles 
0, Q that takes into account the microscopic filament ar- 
chitecture. The Timoshenko theory derives an effective 



FIG. 3: Persistence lengths obtained from the position vari- 
ance. A plot of Eq. ([5]) using the parameters of A = 21 /jm 
and ZJf = 6300 fim obtained by Ref. [J is shown for compari- 
son (solid line). For short microtubules, the data deviate from 
the Timoshenko theory and level to a plateau (dashed). 




2 345 810 15 20 30 

microtubule length L (/im) 

FIG. 4: Drag coefficients ( versus contour length L. For short 
microtubules, the experimental data lie significantly above the 
hydrodynamic estimate £0 of Eq. ([2]) (dashed). The solid line 
shows a fit of Eq. (O taking into account internal friction. 



bending rigidity from stretch and shear terms, which re- 
sult in fully coupled bending for long microtubules and 
an intermediate, shear-dominated length regime where 
rigidity scales with L 2 . In addition to these terms, the 
constituent filaments in a wormlike bundle have a bend- 
ing stiffness that dominates the rigidity of the decoupled 
bundle in the short length regime. 

At the same time, the drag coefficients shown in Fig. [4] 
significantly deviate from the hydrodynamic estimate of 
Eq. (J2|) for short microtubules. The sharp increase in this 
length regime is consistent with the presence of internal 
friction as introduced by Poirier and Marko [2p| . who 
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derived an additional friction term of 
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(I) 



(6) 



to accomodate for dissipation during conformational 
changes or due to liquid flowing through narrow pores 
in a biofilament. In the case of microtubules, both con- 
tributions seem reasonable due to the complex composite 
architecture and the large fluid cylinder inside the fila- 
ment. The parameter e is proportional to an effective 
viscosity summarizing the friction losses inside the fila- 
ment. A fit of the expression 



C - Co + Ci. 



Atttj 



ln(L/d) + 21n(2) - 1/2 



with qi = tt/2 yields a coefficient e — 4.3 ± 0.5 x 
1CP 26 Nsm 2 which is more than one order of magni- 
tude smaller than the numbers of 1.6 x!0~ 24 Nsm 2 and 



6.9 x 10 -25 Nsm 2 obtained by Brangwynne et al. 11 | 
and Janson and Dogterom [lOj, respectively. This dis- 
crepancy suggests that high mode deformations of long 
microtubules may be governed by different dissipative 
mechanisms than first mode deformations of short mi- 
crotubules. 

In this letter, we presented a first analysis of relax- 
ation times and drag coefficients for microtubules in the 
length range most relevant to essential cellular functions 
such as cell division. Our results show that the dynamic 
behavior for this short length regime significantly devi- 
ates from predictions of the wormlike chain model. The 
dynamics seem to be modulated by a length-dependent 
persistence length as well as additional friction contribu- 
tions that may result from internal dissipation. While 
microtubules longer than ~5 /jm display a persistence 
length behavior consistent with the Timoshenko model, 
shorter microtubules show a uniform persistence length 
of 580 ± 100 fim which is consistent with decoupled bend- 
ing of constituent protofilaments as predicted by the the- 
ory of wormlike bundles @, HI • These findings not only 
emphasize the importance of molecular architecture to 
be taken into account for models of semiflexiblc polymer 
mechanics, but also suggest a wealth of possibilities for 
natural modulation of biopolymer properties in cells. 
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